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1 Introduction 

Given a finite subset C C R 2 , calculating the number of triangulations of its 
convex hull Conv(C) using only Euclidean triangles with vertices in C seems 
to be difficult and has attracted some interest, both from an algorithmic 
and a theoretical point of view, see for instance pQ, 0, [3], @], [5], 0; 0; 

The aim of this paper is to describe a class of configurations, convex 
near-gons, for which this problem can sometimes be solved in a satisfactory 
way. Loosely speaking, a convex near-gon is an infinitesimal perturbation 
of a weighted convex polygon, a convex polygon with edges subdivided by 
additional points according to weights. Our main result shows that the 
triangulation polynomial enumerating all triangulations of a convex near- 
polygon is defined in a straightforward way in terms of edge-polynomials 
associated to the "perturbed" edges of such a convex near-gone. These 
polynomials are difficult to compute in general except in a few special cases. 
We present a few algorithms related to them. One of these algorithms is 
a slightly more sophisticated version of an algorithm by Kaibel and Ziegler 
described in |Hj and yields also a general purpose algorithm (unfortunately of 
exponential complexity), for computing arbitrary triangulation polynomials. 
This algorithm, based on a transfer matrix, is fairly simple and it would be 
interesting to compare its performance with existing algorithms, like for 
instance the algorithm of Aichholzer described in pQ . 

Figure^represents a weighted convex triangle and an infinitesimal generic 
perturbation of it. Figure 2 shows an integral realisation isotopic to this in- 
finitesimal perturbation. 

The complete triangulation polynomial of the weighted triangle A de- 
picted on the left side of Figure ^ equals 

901 s 14 + 4825 s 13 + 11734 s 12 + 17130 s 11 + 16710 s 10 + 11466 s 9 
+5670 s 8 + 2034 s 7 + 525 s 6 + 95 s 5 + 11 s 4 + s 3 

showing that there are for instance 17130 triangulations of A using 11 ver- 
tices among the 14 vertices of A. 
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Figure 1: A weighted triangle and an associated generic perturbation 




Figure 2: An integral representant of a convex near-polygon 

This paper contains all details (except for some straightforward calcu- 
lations involving polynomials where computers exceed human possibilities 
by orders of magnitudes) needed for computing the complete triangulation 
polynomial 

194939 s 14 + 338669 s 13 + 263615 s 12 + 119944 s 11 
+34773 s 10 + 6522 s 9 + 748 s 8 + 42 s 7 

of the near-triangle represented in Figure 2. 

In the sequel of this paper, we recall first a few generalities concerning 
finite configurations of R 2 . 

We define then weighted convex polygons and give a formula for their 
complete triangulation polynomial in terms of their edge-weights. 

After introducing convex near-gons, near-edges and edge-polynomials 
we state the main result of this paper expressing the complete triangulation 
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polynomial of a convex near-polygon in terms of the edge-polynomials of its 
edges. 

We list then a few usefull facts and give some data concerning edge- 
polynomials. 

Finally, we describe a transfer-matrix approach for computing edge- 
polynomials (and/or solving enumerative problems concerning triangula- 
tions). 

We close this introduction by describing the following notation, sug- 
gested by the fact that polynomials and formal power series are each other's 
algebraic dual. 

Notations. Given a polynomial p(x) = J2i=o a i x% £ k[a;] and a formal 
power series s(x) = Ya^oPi x% € k[[x]] (over any commutative field or ring 
k) we denote by 

d 

(p, s) x = (p(x), s(x)) x = ^ai Pi 

i=0 

the obvious linear pairing. 

This notation makes sense for polynomials and formal power series over 
several variables and we have 

(P,s t)x, y = ( (P,s)x ,t) y = ( (p,t) y ,s) x 

for a polynomial p G k[x,y] and formal power series s G k[[x]], t G k[[y]]. 
Throughout the rest of this paper, p n denotes always the polynomial 

Pn= E (-!)*( fc J* B - fc eZ[t] 
and C n = { 2 ™)/( n + 1) stands for the n— th Catalan number. 



2 Configurations 

A configuration of the oriented plane R 2 is a finite set C = {P±, . . . , P n } C 
R 2 of n distinct points. Two configurations C = {P\...,P n } and C = 
{P[,...,P^} of n points are (combinatorially) equivalent or isomorphic if 
there exists a bijection <p : C — ► C such that 

sign(det(P j - P h P k - Pi)) = sign(det((^(P,) - ip(P$, ^{P k ) - V (P t ))) 

for all 1 < % < j < k < n where 



sign(x) 



1 if x > 

if x = 
-1 ifx<0 
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Two configurations C = {P\, . . . , P n } and C = {P{, . . . , P^} of n points are 
isotopic if there exists a continuous path t i — ► C(t) G (R 2 )™ of equivalent 
configurations such that C(— 1) = C and C(l) = C. Isotopic configurations 
are equivalent. 

A configuration C C R 2 is generic if three distinct points of C are never 
collinear. 

Given a configuration C, we denote by Conv(C) C R 2 its convex hull. 
An element P G C is extremal if Conv(C \ {.P}) 7^ Conv(C). We denote by 
Extr(C) C C the subset of extremal vertices. They are the vertices of the 
convex polygone Conv(C). 

A triangulation of a finite configuration C is a finite set T = {Ai, . . . , A^,} 
of Euclidean triangles with vertices in C such that Conv(C) = Ua^tA^ and 
Aj n A,- is either empty, or a common vertex, or a common edge of two 
distinct triangles Aj, Aj G T. 

The complete triangulation polynomial p{C) is the polynomial p(C) = 
J2 T k(C) s k G Z[s] whose coefficient Tfc(C) equals the number of triangula- 
tions of Conv(C) using exactly k vertices of C. The number of maximal 
triangulations using all vertices of C is given by the leading coefficient of 
p(C). The lowest non-zero coefficient of p(C) is the (m — 2)— th Catalan 
number C m _2 = C^™-?)/(. m ~ ^) wnere 171 = tt(Extr(C)) is the number of 
vertices of the polygone Conv(C), see for instance Exercice 6.19 of [lflj . It is 
easy to check that the complete triangulation polynomial p(C) depends only 
of the equivalence class of a finite configuration C. 

Remark 2.1 One can also consider the triangulation polynomial defined by 

J2 T koM {C)s^S^ 

counting the number of triangulations using ko vertices and k% edges. The 
number k2 of triangles can then be recovered using the Euler characteristic 
ko — k\ + &2 = 1 of a compact, simply connected triangulated polygonal re- 
gion in R 2 . This more general polynomial yields the same information as 
the complete polynomial considered above except if the boundary d(Conv(C)) 
contains points of C which are not extremal. Most of the results and algo- 
rithms of this paper can easily be modified in order to deal with this more 
general polynomial. For clarity and concision we stick to the simpler version 
p(C) defined above. 

3 Weighted convex polygons 

This section describes weighted convex polygons, a particular set of finite 
planar configurations, and gives formulae for their number of maximal tri- 
angulations and for their complete triangulation polynomials. 
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For I > 3 and natural numbers a\,...,ai > 1 we denote by P{a\ ■ 
a% ■ ■ ■ a{) C R 2 a strictly convex polygon with I (counterclockwise) cyclically 
ordered edges subdivided respectively by (a\ — 1), (02 — 1), . . . , (a/ — 1) ad- 
ditional points. We call such a subdivided polygon a weighted polygon with 
(cyclical) weights a±, . . . , aj. Figure 01 displays two realisations of a weighted 
pentagon P(l • 5 • 2 • 3 • 4). 




(0,5) 



(6,3) 



(0,0) 




(7,2) 



(7,1) 



Figure 3: Two weighted pentagons with edge-weights 1, 5, 2, 3, 4. 

The configurations given by the (ai + . . . + a{) = (61 + . . . + 6/) ver- 
tices of two weighted convex polygons P{a\ ■ ■ ■ a{) and P(b\ ■ ■ -bi) having 
the same number of edges and the same weights up to cyclical permutation, 
are isotopic and their sets of triangulations have the same combinatorics. 
We denote by r^(ai • • • a{) the number of triangulations which use k vertices 
of such a weighted convex polygon P{a\ ■ ■ ■ a{). Such triangulations exist of 
course only if I < k < J2i=i a i- 

For n > 1 we introduce the maximal edge-polynomial p n £ Z [t] by setting 



Pn 




n—k 



eZ[t] 



and the complete edge-polynomial p n € (Z[s]) [t] with indeterminate t over 
the ring Z[s] which is defined as 



p k (t)s" € (ZM) [t] . 



The first few maximal edge-polynomials are 

*5 



Pi 

P2 
P3 
P4 



t 

t 2 -t 

t 3 - 2t 2 

t 4 - 3t 3 + t 2 



P5 
P6 
P7 
P8 



4t 4 + 3t 3 
5t 5 + 6t 4 



V - 6t b + lot 5 - 4r 

t 8 - 7t 7 + 15t 6 - 10t 5 + t 4 



and the first few complete edge-polynomials are 



Pi 
P2 

P3 

Pa 



pi s = t s , 

P2 S 2 + p\ S = (t 2 - t) S 2 + t S , 



p 3 s 3 + 2p 2 s 2 + pi s 



2t z 



P4 S + 3ps S + 3p2 S + pi S 



(t 4 - 3t 6 + 1 2 ) s 4 + 3(t d 



2t z 



■ 3 + 2(t 2 
+ 3(t 2 - 



- t) s 2 + t s , 
t) s 2 +t s . 
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Theorem 3.1 The number r TO ax(oi • &2 • • • o-l) of maximal triangulations of 
the weighted convex polygon P{a\ ■ a<i ■ ■ ■ a{) is given by 

l oo 

<IW*)<E^-2 nt 

3=1 n=2 

where C n = Cn)/( n + 1) * s the n—th Catalan number. 

Corollary 3.2 The complete triangulation polynomial J2k T k( a i ' ' ' o,i)s k of 
a weighted convex polygon P{a\ ■ a% ■ ■ ■ a{) is given by 

l oo 
j=l n=2 

Corollary 3.3 The complete triangulation polynomial (and hence the num- 
ber of maximal triangulations) of a convex weighted polygone P(a± ■ • • oj) de- 
pends only on the multi-set (set having perhaps multiple elements) {oi, . . . , a{\ 
defined by the weights and not on their particular cyclic order. 

Example. The two weighted convex polygons P(l,5, 2,3,4) of Figure El 
have 

(Pi P5 P2 P3 Pi, T,n=2 C n-2 t n ) t 
= ( t ( t 5 _ 4t 4 + 3 ^3 )(t 2 _ t)(t 3 _ 2t 2 ){t A _ 3^ + ^ ^oc^ ^ ^ 

= (t 15 - 10t 14 + 39t 13 - 75t 12 + 74t u - 35i 10 + 6i 9 , £~ =2 C n -2 t n ) t 
= C13 - IOC12 + 39Cu - 75Ci + 74C 9 - 35C 8 + 6C 7 
= 7429000 - 10 • 208012 + 39 • 58786 - 75 • 16796 

+74 • 4863 - 35 • 1430 + 6-429 
= 8046 

maximal triangulations. 

Their complete triangulation polynomial J2k T k{P{^ • 5 • 2 • 3 • 4))s fe equals 

(Pi P~5 P~2 P3 P~4, T,n=2 C n - 2 t n ) t 

= 8046 s 15 + 37250 s 14 + 77467 s 13 + 95364 s 12 + 77048 s 11 

+42776 s 10 + 16584 s 9 + 4460 s 8 + 805 s 7 + 90 s 6 + 5 s 5 . 




Figure 4: A maximal triangulation with 7 (shaded) up-triangles of the digon 
P(9-7). 
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Remark 3.4 Th,eorem, \3.1\ and Corollary \S. Si remain valid for weighted "di- 
gons", weighted "convex" polygons P{a\ ■ 02) having only two edges. Their 
"sides" or "edges" are in this case not straight lines but piecewise linear 
paths bending sligthly outwards in order to enclose a non-empty interior if 
a\ 02 > 2 (see Figure 4 f or an example). The weighted digon P(l, 1) has 
by convention one triangulation. P(0, a) = P(a, 0) has 1 triangulation for 
a = and none otherwise. In all the remaining cases, only "triangulations" 
with all triangles meeting both edges are allowed. The number of such tri- 
angulations of P{a\ ■ 02) equals ( ai a^-2~ 4 ) : Indeed, every triangle except the 
first and the last one is either up (with two vertices on the upper edge) or 
down (with two vertices on the lower edge) and there are exactly a± — 2 such 
triangles which are up, cf. Figure 4 for an example. 



3.1 Proofs for weighted convex polygons 

Let -P(a" 1 • a^ 2 ■ ■ ■ af) denote the weighted convex polygon having a\ edges 
of weight a\ followed by a 2 edges of weight 02 etc. We denote by r m ax(«i 1 ■ 
a 2 2 ' ' ' a T) ^ ne number of maximal triangulations of P^a® 1 ■ ■ ■ af l ). 

The main ingredient for proving Theorem l3.1l is the following illustration 
of the Inclusion- Exclusion Principle (cf. Chapter 2 of of 

Proposition 3.5 (Inclusion-Exclusion Principle.) We have 

Lai/2J (a x -k\ 
rmax(ai ■ a 2 ■ ■ ■ ai) = X^ -1 )^ k J T max(i ai k -a 2 ---a l ). 

Proof. Enumerate the sides of P(l ai • 02 ■ ■ ■ a;) cyclically such that the first 
a± sides are of weight 1 and correspond to the initial factor l ai of the product 
l ai ■ 02 • • • a>i- Let T be a maximal triangulation of the weighted convex 
polygon P(l ai ■ 02 • • • ai). An initial boundary triangle of T is a triangle of 
T having two edges contained among the first a% sides (having weight 1) of 
the convex weighted polygon P(l ai ■ 02 ■ ■ ■ a{). A subset {5\, . . . , 5k}, < 
k < L a i/2j °f k initial boundary triangles in T is called a A ^— decoration 
of T. Starving to death the k triangles 5\,...,5k of the A^— decorated 
maximal triangulation (T, {5\, . . . , 5k}), we get a maximal triangulation T of 
P(l ai_fc • 02 ■ ■ ■ ai) together with a graveyard {ei, . . . , e^} of k marked edges 
in memory of the deceased triangles. We call the pair (T, {ei, . . . , e^}) a 
maximal E^— decorated triangulation of P(l ai ~ k -02 ■ ■ ■ ai). Figure 5 displays 
a A2— decorated maximal triangulation of P(l 7 -3-2) and the corresponding 
E2— decorated maximal triangulation of P(l 5 -3-2). 

Obviously, A^— decorated maximal triangulations of P(l ai -02 ■ ■ ■ ai) and 
E^— decorated maximal triangulations of p(l ai_fc • a 2 • • • ai) are in bijection 
and since any maximal triangulation of P(l ai ~ k -02 ■ ■ ■ ai) has (^\ k ) different 
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Figure 5: Decorated maximal triangulations of P(l 7 -3-2) and P(l 5 -3-2) 
Efr— decorations, there exist exactly 



fax - k\ 

I * ) 



rmax(l ai • «2 • • • oj) 

maximal decorated triangulations of P(l ai ■ a>2 ■ ■ ■ ai). 

Consider now a maximal triangulation T of P(l ai ■ a>2 ■ ■ ■ a/) having ex- 
actly s initial boundary triangles. The triangulation T gives rise to (f) dif- 
ferent A/%— decorated maximal triangulations of P(l ai ■ a,2 ■ ■ ■ a/) and yields 
a contribution of 




-1)' 



1 if s = 
otherwise 



to the alternating sum 



k 




rmax(l ai k a 2 ■ ■ ■ a/) 

counting decorated maximal triangulations of P(l ai -a^ ■ ■ - ai) with sign 
(— l) k . This alternating sum counts thus exactly the number of maximal 
triangulations of P(l ai • 02 • ■ ■ a{) without initial boundary triangles. Such 
triangulations are in bijection with triangulations of P(a\ -ai ■ ■ ■ a\): straight 
out the first a\ edges of weight 1 in P(l ai • 02 • • • aj) into a unique first edge 
of weight a\ of P{a\ ■ ■ ■ a{). Figure l3~T1 illustrates this: Its left side displays a 
maximal triangulation without initial boundary triangles of P(l 5 -3-2) and 
its right side shows the corresponding straightened maximal triangulation 
ofP(5-3-2). □ 
Proof of Theorem 13.11 Given a polynomial q(t) = J2i=o a i t % we se t 

d 

T(q(t); a 2 - ■ ■ ai) = ^ai r m ax(l 4 • a 2 • • • a/) . 
i=0 

Since the polynomial 

L»i/2J 

Pa '' 



(t)= E(-i) fc ( k y- 
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Figure 6: Straightening a maximal triangulation without initial boundary 
triangles 



mimicks the formula of Proposition 13.51 we get by iterated application of 
Proposition 13.51 

i 

rmax(ai ■a 2 ---ai) = r(J^[p ai (i); 0) = r max (l' J ) 

i=l j 

where Y^jljV = Ui=iPaM- 

Since P(l n ) is a convex polygon with n edges (of weight 1) and n ex- 
tremal vertices, the number Tmax(l n ) of maximal triangulations of P(l n ) is 
given by the (n — 2)— th Catalan number C n _2 = ( 2 fc^)/fa ~ -0> ^ or 
instance Exercice 6.19 of ^Sl- Hence the result. □ 
Proof of Corollary 13.21 Non-maximal triangulations of Piflx ■■■ai) are 
in bijection with maximal triangulations of P{b\ ■••&/) where 1 < b{ < a,; 
and where the weighted polygons P{b\ •••&;) are realised by choosing in all 
(b'Zi) possible ways subsets of (bi — 1) points among the (a^ — 1) interior 
points of the i— th edge having weight en of P(a± ■ ■ ■ a{). 

The complete triangulation polynomial of P{a\ ■ ■ ■ ai) is thus given by 

l<bi<a, \j=l \" 1 V j j=l n=2 
1 aj ( — 1 \ 00 

j=l 6 J= 1 \°J V n=2 
I 00 

j=l n=2 

which proves the result. □ 
Corollary 13.31 is obvious. 



4 Near-gons 

An n— near-edge E is a sequence of (n + 1)— points 



(P ,...,P n ) = (\ 



X \ / XX \ I X n -l \ I x n 

vo r\vi r'"' \ yn-i r\vn 
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in R 2 such that xo < x\ < . . . < x n _i < x n . We consider the complete 
order Pq < Pi < P2 < . . . , P n on E and call Pq, respectively P n , the initial, 
respectively final, vertex of the near-edge E. We denote a near-edge E either 
by the sequence E = (Pq, . . . , P n ) of its points or by the 2 X (n + 1) matrix 



E 



satisfying xq < x\ . . . < x n whose columns contain the coordinates 



of Pi. We endow E with the total order Pq < Pi < . . . < P n . 

Two n— near-edges E = {Pq, . . . ,P n } and E' = {Pq, . . . , P^} are equiva- 
lent if 

sign (det(P j - Pi, P k - Pi)) = sign (det(Pj - P[, P' k - Pi)) 
where 0<i<j<k<n and where sign(x) <G {±1, 0} is defined by 



sign(x) 



1 if x > 

if x = 
-1 ifx<0 



A continuous path t 1 — ► E(t), t G [—1, 1] of equivalent near-edges is an 
isotopy. Two near-edges E{— 1) and E(l) joined by an isotopy are isotopic. 
Isotopic near-edges are of course always equivalent. 

A near-edge is generic it the underlying set of points is a generic config- 
uration of R 2 , i.e. is without three collinear points. 

Given a near-edge E = (Pi = ^ ^* ^ ) i=0 n and e E R, we denote by 

E e the near-edge 

E " - (p - (e) - ( T Vi ) - ( \ \ ) • 

For e > 0, the near-edges E and E e are isotopic and thus equivalent. 

Let P be a convex polygon with I extremal vertices {Vo, V\, V2, ■ ■ ■ , VJ_i, V\ = 
Vq} appearing in counterclockwise order around its boundary dP. 

Given a sequence E^ , . . . , E\ l where Ei = (Pj,o, Pi,i, ■ ■ ■ , P%,m) 1S a n i~ near- 
edges and £j£Ra real number, we denote by 

P{E?-E?...E?) 

the unique configuration obtained by gluing the n^— near-edge E\ l , rescaled 
by a suitable orientation-preserving similitude, along the oriented edge of P 
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which starts at V%-\ and ends at Vi- More precisely, the gluing similitude ipi 
is the unique orientation-preserving similitude of R 2 such that 

The configuration P{E^ ■ ■ ■ ■ Ep) is now the set of points l4 =1 (/?j(£^) C 
R 2 . 

We have the following result which we state without proof. 

Proposition 4.1 (i) The configurations PiEl 1 ■ ■ ■ Ef l ) are isotopic for < 
6{ small enough. 

(ii) Given a second convex polytope P' having I vertices Vq , V[, . . . , V{ = 
Vq in counterclockwise order, the configurations 

PiEl 1 ■ ■ ■ El 1 ) and P'(E? ■ ■ ■ E e L l ) 

associated to P and P' are isotopic for < 6j small enough. 

(Hi) Given pairs (E^E^) of equivalent near-edges, the configurations 

P(E?---E?) and P{E[ ei ■ ■ ■ E[ ei ) 

are equivalent for < e, small enough. 

We call the equivalence class of the unique configuration described by 
Proposition 14.11 fwith < small enough) the convex l—near-gon or near- 
gon with (cyclically oriented) near-edges Ei, . . . , E[. We denote by P{E\ • ■ ■ E\) 
any configuration representing the equivalence class of this near-gon. 

Theorem 4.2 For every near-edge E there exist polynomials 

p(E) € Z[t] andp(E) G Z[s,t] 

such that the convex near-polygon P(E\ ■ ■ ■ E{) has 

l oo 

{Y\p{E 3 )^c n -2 n t 

j=l n=2 

maximal triangulations and complete triangulation polynomial 

I oo 
j=l n=2 

where C n = ( 2 ^ 1 ) /(n + 1) are the Catalan numbers. 

Corollary 4.3 The number of maximal triangulations and the complete tri- 
angulation polynomial of a convex near-polygon P{E\, . . . ,Ei) depend only 
on the multiset {E\, . . . ,Ei} of (equivalence classes of) its near-edges. 
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Example. The configuration of Figure 2 is equivalent to the convex near- 
gon P(E a ■ E b ■ E c ) where 




are near-edges whose complete edge-polynomials 

p(E a ) = (14p 3 + 7p4+P5)s 5 + (10p 2 + 7p3 + 2p 4 )s 4 

+(2pi + 2p 2 +p 3 )s 3 , 
p(E b ) = (5p 3 +p 4 )s 4 + 2(2p 2 +P3)s 3 + (pi+p 2 )s 2 
p(E c ) = (10p 3 + 7p 4 + 2p 5 )s 5 + (3p 2 + 13p 3 + 4p 4 )s 4 

+3(2p 2 +p 3 )s 3 + (pi+p 2 )s 2 

(where p n = El={? J (-l) fc ( n fc fc )^~ fc G Z M as usual ) wil1 be computed in 
the sequel. The complete triangulation polynomial p(E a ■ E b ■ E c ) of this 
configuration is given by 



(p(E a )p(E b )p(E c ),j:n=2Cn-2 t n ) t 

= 194939 s 14 + 338669 s 13 + 263615 s 12 + 119944 s 
+34773 s 10 + 6522 s 9 + 748 s s + 42 s 7 . 



4.1 Edge-polynomials 

Given an n— near-edge 

*-«..-.*.>-(;: : ::: L") 

we define the lower convex hull d-E as the piece- wise linear path joining Pq 
to P n formed by the "lower" edges of <9Conv(.E). The set V-(E) of lower 
extremal vertices of E is defined as the increasing subsequence ~Extr(E)Pid~E 
of all extremal vertices of Conv(^) situated in the closed "lower" halfplane 
containing Pq and P n in its boundary and containing all points (xo,y) with 

y < yo- 

A roof of E is a strictly increasing sub-sequence 

R = ( P jo = -fb, Pji , ■ ■ ■ , Pj k -i ) -Pjfc = Pn) C E 

starting at Po and ending at P n of E 1 . If a roof contains k + 1 elements 
we call length(i?) = fc its length. The skyline Skyline(i?) of a roof i? is the 
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piecewise- linear path from Pq to P n defined by joining consecutive elements 
of R using straight segments. 

A roof R is covering if it "shelters" every element of E: for every point 

Pi = ( ^ ) G E W6 haV6 either p i e R or Vi > Vi for X ~. j G Skyline(P). 
The maximal edge-polynomial p(E) of the near-edge E 1 is defined as 

P( E ) = E T ma x(E,R) Pi engt h(R) G Z M 

H covering roof of E 

where T m ax(P) P) denotes the number of maximal triangulations of the 
(generally non-convex) compact polygonal region delimited by the two piece- 
wise linear paths Skyline(P) and d-E. 

A k— sub-edge of E is an increasing subsequence E'cEofk+l<n + l 
elements in E such that V-(E') = V-(E). In particular, any sub-edge E' of 
E has also initial vertex Pq and final vertex P n . 



Figure 7: An 8— near-edge and a covering roof of weight 4 of a 6— sub-edge 
Example. The left side of Figure 7 displays the 8— near edge 

/ 123 4 5 6 7 8 \ 



E=(P ,...,P 8 



0-111-2-3 -2 -1 



We have V-(E) = (P ,Pi,P 5 ,P 8 ) and E has 2 5 sub-edges obtained by re- 
moving any subset of vertices among {P%, P3, P4, Pg, P7} from P. The right 
side of Figure 7 displays the covering roof R = (Po, P3, P5, P7, Ps) of the 
sub-edge P' = {P , P u P 3 , P 4 , P 5 , P 7 , P 8 ) C P. 

The complete edge-polynomial p(E) of an n— near-edge E = (Po, . . . , P„) 
is defined as 

n 

p( E ) = E E p(^')« m 

m=1 E'cE m-sub-edge 

where p(E') denotes the maximal edge-polynomial of the sub-edge E' C E. 

Example. We compute the complete edge-polynomial p(E a ) of the near- 
edge 

P o = (P ,...,P 5 )= ^ 1 J ^ -1 ) 
involved in the near-gon P(E a ■ E^ ■ E c ) of Figure 2. 
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Figure 8: All covering roofs of sub-edges involved in p(E a ) 



Figure 8 contains all roofs of the four possible sub-edges of E a obtained 
by removing any subset of points in {Pi,P 3 }. Each possible sub-edge has 
four different covering roofs whose contributions to p(E a ) are given by 



sub-edge 










(Po, Pi, P2, P3, Pii -P5) : 


14p 3 s 5 


2p 4 s 5 


5p4S b 


p 5 s 5 


(Po, P2, P3, Pi, P5) '■ 


5p 2 s 4 


2p 3 s 4 


2p 3 s 4 


4 


(Po, Pi, P2, P4, P5) '■ 


5p 2 s 4 


4 




4 


(P ,P2,P4,P 5 ) ■ 


2p 1 s' i 


p 2 s :i 


P2S 3 


p 3 s 3 



They sum up to the complete edge-polynomial 

p(E a ) = (Up 3 + 7p 4 + p 5 )s 5 + (10p2 + 7p 3 + 2p A )s 4 + (2p! + 2p 2 + p 3 )s 3 
of £«. 

4.2 Proof of Theorem I4T21 

Consider a triangulation T of a convex near-polygon P{E\ ■ ■ ■ E{) with com- 
plete triangulation polynomial 

k 

A triangle A G T with vertices {a, 6, c} is of edge-type if there exists j such 
that {a, b, c} C J5j. Otherwise, the triangle A is interior. If A with vertices 
{a, b, c} is of edge-type, we write A E Ej if {a, 6, c} C Ej. 

The definition of a convex near-polygon implies that the set of all tri- 
angles A € Ej n T determines a covering roof Rj of a suitable fcj— sub-edge 
Ej of Sj. More precisely, define Ej as the set of all (kj + 1) points of Ej 
which appear as vertices in T. The set of edge-type triangles in Ej tri- 
angulates now a (generally non-convex) polygonal region enclosed by two 
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(non-crossing) paths joining the endpoints of Ej. These paths are the lower 
convex hull d—Ej = d-Ej and the skyline Skyline(-Rj) of a unique covering 
roof Rj C Ej having length length(i?j). The interior triangles of T define 
(after straightening) a unique maximal triangulation of the weighted convex 
polygon P(length(i?i) • • • length (Ri)). The left side of Figure 9 displays a 
maximal triangulation of a near-gon having four near-edges (of length 2,3,1 
and 4). The inner triangles of the displayed triangulation yield a maximal 
triangulation of the convex weighted polygon P(2 • 2 • 1 • 3) depicted on the 
right side of Figure 9. 




Figure 9: Edge- type and interior triangles of a triangulation 
There are exactly 

Tmax (length^) • • • length^)) 1Q 

Tmax(-Ej, Rj) 

3=1 

such triangulations contributing to the coefficient s i k i+ k z+---+ k i) Q f p for fixed 
covering roofs Rj of fixed kj— sub-edges E'a C Ej. Applying Theorem 13.11 
and summing over all covering roofs in sub-edges yields the result. □ 

5 Properties of edge-polynomials 

A near edge E = (Pq, . . . , P n ) factorises into near-edges Ei,E 2 if there exists 
a lower extremal vertex P^ £ E n d-E such that 

El = (Po, Pi, ■ ■ ■ , Pk~l,Pk), E 2 = {Pk, Pk+l, • • , Pn-l,Pn) 

and all points of E \ Ei lie strictly above any line defined by two distinct 
points of Ei for i = 1,2. We write E = Ei ■ E% if the near-edge E factorises 
with first factor Ei and second factor E 2 . A near-edge is prime if it has 
no non-trivial factorisation. It is easy to show that every near-edge has a 
unique factorisation into prime near-edges. 

Proposition 5.1 Given a factorisation E = E% - E 2 of a near-edge E we 
have 

p{E)=p{Ei) p{E 2 ) andp(E)=p(Ei) p(E 2 ) . 
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Proof. Since the near-polygons P(l k ■ E) and P(l k ■ E\ ■ E2) are equivalent 
for all k = 2, 3, . . . we have 



(t k p(E), e c n _ 2 n t = (t k pis,) p(e 2 ), jr c n _ 2 r) t . 

n=2 n=2 

This implies the result since det((Ci + j + k) 0<i j <n ) > for all k > and 
n > 1 (this follows for instance easily from Exercice 6.26.b in 1=1 
Since the so-called A;— shifted Hankel matrices 

Hi j = C i+j+k , < i,j < n 

are all non-singular, the coefficients 7 Q ,...,7d (with a > 2) of a polyno- 
mial = X)f= a 7i ^ can be obtained from the (d — a + 1) "linear form 
evaluations" 

00 

{t k q(t),J2Cn-2 t n ) u k = 0...,d-a. 

n=2 

This turns any method or algorithm for enumerating maximal (or all) 
triangulations of configurations into a method or algorithm for computing 
maximal (or complete) edge-polynomials. 

Example. Using the online computing service [S] of O. Aichholzer, we 
compute the maximal triangulation polynomial p(E c ) (which is of of course 
also the coefficient of s 5 in the complete triangulation polynomial p{E c )) of 
the near-edge 

' 1 2 3 4 5 \ 
2 1 -1 1 J 



involved in the near-gon of Figure 2. 

Since the set d + (E c ) of all boundary vertices on the upper boundary of 
Conv(i? c ) consists of the 3 segments P0P1, P1P4 and P4P5, we have 

p(E c ) = ap 3 + (3p4 + 7p 5 

where a, (5 and 7 are unknown natural integers. 

We identify now the near-edge E c = (Pq, . . . , P5) with its underlying set 



Define 



and check that the configurations 

^ c U{^i}, E c U{Ai,A 2 }, E C U{A 1 ,A 2 ,A 3 } 
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are near-polygons with near-edge factorisation 

E c ■ E f , E c ■ E f , E c ■ E\ 

where E\ denotes the unique 1— edge with maximal edge polynomial p{E\) = 
pi = t. Denoting by T max (C) the number of maximal triangulations of a 
finite configuration C and using for instance the on-line computing service 
of Aichholzer mentionned above we get the linear system of equations 

19 = W4U{il}) = (p(E c )t 2 ,Y J n=2Cn-2t n ) t 

87 = r ma x(^U{ii,i 2 }) = <p(E c ) i 3 ,E£°= 2 C n -2 t n ) t 
175 = T max (E c U{A 1 ,A 2 ,A 3 }) = (p(E c )t\J2n=2C n -2t n )t 

Substituting p n = El=o J (-l) fc ( n fe C n = ( 2 n n )/(n + l) and solving for 

the unknowns a, /3, 7 we get 

p(E c ) = 10p 3 + 7 P4 + 2p 5 . 

Remark. Since the computations of the coefficient a = T max (-E'c) = 10 of 
lowest order (in t) and of the leading coefficient 

7 = T ma *({(o J,( I ).(?).( _l)» ' 

■W{(_?).(j).(5)»-2-l-2 

are easy, we could have simplified the above computation by a considerable 

amount of work. 

For completeness we mention also the follwing obvious fact: 

Call two n— near-edges E = (Pq, . . . , P n ) and E' = (Pg, . . . ,p' n ) vertical 

mirrors if E' given by E' = (P n , P n —i, ■ ■ ■ > Pi: Po) where Pj = I Xl \ 

\ Vi J 

is the Euclidean reflection of Pi = ^ X% ^ with respect to the vertical line 
x = 0. 

The following result is obvious: 
Proposition 5.2 If a pair of near- edges (E,E') are vertical mirrors, then 

p{E) = p(E') . 

5.1 Complete edge-polynomials for small near-edges 

This subsection contains representants and edge-polynomials for all 1—, 2— 
and 3— near-edges. 
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5.1.1 1— near-edges 

The unique 1— near-edge can be represented by E± = ^ ^ ^ ^ . It is generic 
and prime and has complete edge-polynomial p{E\) = p\s = s t. 

5.1.2 2— near-edges 

There are two generic 2— near-edges, represented by 

E 21 = ( ° 1 2 ) E 22 =( ° 1 2 ) 
2A \0 1 J ' 2 ' 2 ^ -1 o ) ■ 

E 2 ,i is prime while E 2>2 = E\ ■ E\. They have complete edge-polynomials 

p(£ 2 ,i) = P2S 2 + Pis, p(E 2;2 ) = (p 2 + Pl )s 2 = piExf . 

Moreover, there is also a unique non-generic 2— near-edge represented for 
instance by 

^ 3= (o o o) 

with complete edge-polynomial given by 

p{E 2 , 3 ) = P 2 s 2 + pis = p(E 2)1 ) . 

5.1.3 3— near-edges 

There are (up to equivalence) eight generic 3— near edges represented by 

, 1 2 3 \ /0 1 2 3\ 

3,1 ~ 1 1 3 J ' 3 ' 2 ~ \ 1 1 ) ' 

, 1 2 3 \ ( 1 2 3 \ 

1 ' 1 ! 3 1 ) ' 3 ' 4 ~ I 1 -1 ) ' 

M _ , 1 2 3 \ / 1 23 

() 1 1 J ' 3 ' 6 ~ I -3 -1 

1 2 3 \ _ / 1 23 

3 ' 7 ~ 1 -1 -1 J ' 3 ' 8 ~ ^ -1 -3 

The first five are prime. The last three have factorisations 

7T1 TP TP TP Tp3 TP TP TP 

^3,6 — ^1-^2,1 ,-^3,7 — t^li ^3,8 — ^2,1-^1 • 

The pairs {E 3j1 ,E 3j3 }, {E 3A ,E 3:5 }, {E 3fi ,E 3j8 } are vertical mirrors. The 
prime near-edges have complete polynomials 

p(E 3 ,i) = p{E 3)3 ) = (p2 + p 3 )s 3 + 2p 2 s 2 + Pis , 

p{E 3>2 ) = 2p 3 s 3 + 2p 2 s 2 +pxs , 

V{E 3A ) = p{E 3 , 5 ) = (2p 2 +p 3 )s 3 +pis 2 . 
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There are moreover nine more 3— near-edges which are not generic. They 
are represented for instance by 







1 


2 


CO 








1 








1 


2 


3 





1 











1 


2 


3 





2 


1 








1 


2 


3 


















1 


2 3 





-1 


-2 





1 


2 3^ 








-io; 





1 


2 3} 





-1 


j 





1 


2 3 





-2 


-1 



Es ' 9 - { 1 2 ) ' Es > 10 
F - ( ° 1 2 3 ) F - 

£3,11 - I o o i o ) ' E3 > 12 ~ 

F 0123 F — 

#3,13 " ( o x o ) ' E3 > lA ~ 

F - ( ° 1 2 3 1 F 

^3,15 — I n o 1 n I > ^3,16 
-£■3, 17 = 

The near-edges 

TP TP IT TP TP TP 

-^3,10 — ^2,3 -C'lj -^3,16 — -C'l ^2,3 

factorise. The remaining prime near-edges have complete edge-polynomials 

p(E 3>9 ) = p(E 3)15 ) = p 3 s 3 + 2p 2 s 2 +pis, 
P(E3,ll) = p(E 3) i 3 ) = (p 2 +P3)s 3 + 2p 2 s 2 +pis, 

P{E 3 ,1 2 ) = P( E 3,14) = (P2+P3)s 3 +P2S 2 , 

V{E 3 ,n) = P 3 = p 3 s 3 + 2p 2 s 2 + pis . 



5.2 Enumeration of generic near-edges 

This subsection is a digression sketching briefly a bijection between near- 
edges (up to equivalence) and orbits of suitable extremal vertices of config- 
urations under automorphisms. 

Let P € Extr(C) be an extremal vertex of a configuration C = {Pi, . . . , P n +2] C 
R 2 consisting of (n + 2) points. Suppose that P is not collinear with two 
distinct elements of C \ {P}- Let P-,P+ € Extr(C) be the immediate pre- 
decessor and successor (in counterclockwise order) of P on the boundary 
dConv(C). 

A suitable projective transformation (sending a line outside Conv(C) 
which misses P only nearly to the line at infinity and sending P close to 
the ideal point (0, +00)) transforms C \ {P} into an n— near-edge E with 
initial vertex P + and final vertex P_. This near-edge E is well-defined, up 
to equivalence, and every n— near-edge can be obtained in this way. 

In particular, generic n— near-edges are in bijection with orbits under au- 
tomorphisms of extremal vertices in generic configurations of (re + 2) points. 
Figure 10 illustrates this for n = 3. It shows all three generic configurations 
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Figure 10: Generic configurations of 5 points and the associated generic 
3— near-edges 

of five points with extremal vertices labelled by the corresponding generic 
3— near-edges (using the notations of the previous subsection). 

Call two generic n— near-edges related if they are associated to two ex- 
tremal vertices of the same generic configuration of (n + 2) points. Classes 
of related generic n— near-edges are in bijection with equivalence classes of 
generic configurations having (n + 2) points. 

We say that two configurations E = (Pq, ... ,P n ) and E' = (Pq, . . . , P^) 

are related by a horizontal reflection if P[ = ( X * ^ is the reflection with 

respect to the horizontal line y = of Pi = 

Equivalence classes of related generic n— near-edges, up to horizontal 
reflections, are in bijection with generic projective configurations (up to 
equivalence) of (n + 2) points in the projective plane RP 2 . 

6 An algorithm for convex near-edges 

An n— near-edge E = {Po, . . . , P n } is convex if its points Po,...,P n are 
all extremal, i.e. vertices of a convex polygon with (n + 1) edges. This 
section describes an algorithm for computing maximal and complete edge- 
polynomials of convex near-edges. 

A convex n— near-edge can be represented by the sequence of points 

"-(;> 

where e±, . . . , e ra _i G {±1}- There are 2™ _1 equivalence classes, encoded by 
(ei, . . . , e n _i) G {±l} n_1 , of convex n— near-edges. The convex near-edge 
with 6i = — 1, i = 1, . . . , (n — 1) has the factorisation Ef . All others are 
prime. 

The following algorithm of polynomial complexity computes the com- 
plete (respectively maximal) edge-polynomial of a convex near-edge. It is 
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based on the fact that all connected components of the open polygonal re- 
gion delimited by a covering roof R and d-(E) are open convex polygons 
if E is a convex near-edge. The exponent of u keeps track of the length of 
a "covering roof under construction" (Pq, P^ , . . . , Pj . , ?). The exponent of 
w keeps track of the number of vertices of the rightmost polygonal region 
which has not yet closed up. 
Algorithm 

Input: n £ 1, 2, ... , and (ei, e 2 , . . . , e n _i) G {±l} n_1 . 
Set 

{s ■ u (complete edge-polynomial) 

u (maximal edge-polynomial) 

For i from 1 to n — 1 do { 
If £i = 1 replace 

{r • (s ■ u ■ w + 1) (complete edge-polynomial) 
r ■ u ■ w (maximal edge-polynomial) 

If ei = — 1 replace 



r ■ s ■ w + ( r, { ^ C n u> n ) ) w • s • i* (complete edge-polynomial) 



vra=0 



oo 



r ■ w + (r, { ^2 C n w n ) w u (maximal edge-polynomial) 

\n=o / 



} 

Result: 



Vn=0 




(where p n = Efc(-l) fc ( n fc fc )* n_fc )- 
Example. The first two near-edges 



1 2 3 4 5 
1-11-10 



and 

' 1 2 3 4 
1-110 



involved in Figure 2 are convex. The near-edge E a corresponds to n = 5 
and the sequence (1, —1, 1, —1). The near-edge E^ corresponds to n = 4 and 
(1,-1,1). 

The above algorithm for E a (with n = 5 and (ei, . . . , 64) = (1, —1, 1, —1)) 
yields 
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i = 


Ci = 


r = 






su 


1 


1 


s 2 u 2 w + su 


2 


-1 


s 3 w 2 it; 2 + s 2 uw + s 3 u 3 + s 2 u 2 


3 


1 


s u 3 w 3 + s 3 u 2 w 2 + s 4 n 4 w + s 3 u 3 w+ 






+s 3 u 2 w 2 + s 2 «w + s 3 u 3 + s 2 u 2 


4 


-1 


s 5 w 3 it; 4 + 2s 4 u 2 u? 3 + s 5 u 4 w 2 + s 4 u 3 w 2 + s 3 -uu> 2 + s 4 ii 3 tt; 






+s 3 u 2 w + 5s 5 u 4 + 4s 4 u 3 + s 5 u 5 + 2s 4 u 4 + s 3 u 2 + s 3 u 3 



Result: (14p 3 + 7p A + p 5 )s 5 + (10p 2 + 7p 3 + 2p 4 )s 4 + (2 Pl + 2p 2 + P3)s 3 

which is the complete edge-polynomial p(E a ) of £7 a . 
Computing 

<r, >«,» 

with 

r = s 4 u 3 tt; 3 + 2s 3 u 2 w 2 + s 4 u 4 w + s 3 u 3 it; + s 2 uu> + s 3 u 3 + s 2 u 2 
corresponding to i = 3 we get the complete edge-polynomial 

p(E b ) = (5p 3 + Pa)s a + (4p 2 + 2p 3 )s 3 + (p! + p 2 )s 2 
of the convex near-edge Ef, encoded by n = 4 and (1, —1, 1). 

7 Roof-sequences and a transfer matrix 

In this section, we define the set of roof-sequences of a finite configuration 
(or near-edge). Such sequences are in bijection with triangulations. Using 
this bijection, we introduce a transfer matrix for computing the complete 
triangulation polynomial (or complete edge-polynomial) of a configuration 
(or near-edge). 

We endow R 2 with the total order relation defined by 




if either x < x' or x = x' and y > y' . This order is obtained by reflecting 
(orthogonaly) the plane endowed with the familiar lexicographic order with 
respect to the line y = 0. 

Consider a finite configuration C = {Pq, . . . ,P n } C R 2 consisting of 
(n + 1) totally ordered points Pq < P± < . . . < P n for the above order. We 
identify C with the strictly increasing sequence (Po, . . . ,P n ) of its (n + 1) 
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points. The initial point Pq and the final point P n of C are extremal vertices 
of Conv(C) and are joined by two (generally distinct) piecewise-linear paths 
<9_l_(C) and <9_(C), called the upper and lower convex hull. The path d+(C) 
(respectively <9_(C)) is contained in the boundary 9Conv(C) and links Pq to 
P n by turning clockwise (respectively counterclockwise) around the convex 
hull Conv(C) of C. 

A roof R of C having length w = length (R) is an increasing subsequence 
R = (P io = P ,P h , . . . , P iw _ 1 , P iw = P n ) containing (w + 1) elements. 

The skyline of a (decorated) roof R = (Pq, ... ,P n ) is the piece- wise lin- 
ear path from Pq to P n obtained by joining consecutive points of R using 
rectilinear segments. Skylines and roofs are in bijection for generic configu- 
rations. If C is non-generic, several roofs might share a common skyline. 




r 10 

Figure 11: A triangulated configuration of (10 + 1) points 
Example. The configuration 



C — {Pq, • • • , P n } 



0000122233 



3210442132 -1 / 

represented in Figure 11 has four roofs with skyline cLC: 

(■Fb, P3, Pio) , (Po, Pi, P3, Pio) , 

(Po, P2, P3, Pio) , (Po, Pi, P2, P3, Pio) 

and two roofs with skyline d + C: 

{Po,Pi,Pb,P&,Plo) {Po,Pi,Pb,P&,P$,Plo) ■ 

A decorated roof is a roof 

P = {Pio = Pq, Pi\ , ■ ■ ■ , P'iw-i , P'im = Pn) , 

decorated by one of the segments joining two consecutive points Pi d ,Pi d+1 
of R. We denote such a decorated roof by 



R — {Pq, ■ ■ ■ , P%d-i , P'id ' Pid+l ' Pid+2 , ■ ■ ■ , Piw — Pn 
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and call R initially decorated if d = 0, i.e. if the decoration starts at the 
initial point Po 

Three non-collinear points A < B < C £ C are vertices of a triangle 
0/ type V if B is "above" the line containing A and C. The non-collinear 
points A, B, C are vertices of a triangle of type A otherwise. In other terms, 
a triangle with (necessarily non-collinear) vertices A < B < C is of type A 
if det(C - A, B - A) > and of type V if det(C - A, B - A) < 0. 

We call a triangle with non-collinear vertices A,B,C minimal if {A, B,C} = 
C n Conv(^4, B,C). In particular, a triangulation is maximal if and only if 
all its triangles are minimal. 

A decorated roof 

P = {Rio = Ro, ■ ■ ■ i Rik-i J R'ik 1 Rik+i ' Rh+2 ' Rik+3 1 ■ ■ ■ 1 Riw = Rn) 

of length w is a A— successor of P if R is a decorated roof of length w — 1 
given by 



^ — (-Pio — -fb) ■ ■ ■ j -fid j -f«d+i ) • • • > Rik-i ' ^«fc ' ^»fc+2 ' Rik+3 ' • • • ' ^««> — 

satisfying Pj d < Pj fe and the points {Pj fe , Pi fc+1 , Pi k+2 } C R' are vertices of a 
A— triangle. 

Otherwise stated, R' is obtained from P by gluing a A— triangle with 
decorated left-upper edge onto an edge of R which is not situated to the left 
of the decorated edge in P. 

Similarly, a decorated roof 



P — (Pi — Po> ■ • • j Rik-i 1 Rik ' Rik+2 ' Rik+3 ' • • • ' — 

of length w — 1 is a successor of P if P is the decorated roof of length w 
given by 



P — (Pio — PO) ■ ■ ■ j Pi d t Pi d +i 1 ■ ■ ■ 1 Rik ' Rik+l ' Rik+2 ' ^«fe+3 ' • • • ' Riw ~ Rn) 

satisfying P id+1 < P ik+2 and the points {P fe , P ik+1 , P ik+2 } C P are vertices 
of a y— triangle. 

Otherwise stated, R' is obtained from P by gluing a V— triangle with 
decorated upper edge onto two consecutive edges of increasing slope of P and 
the decorated edge of P is not situated to the right of these two consecutive 
edges. 

Example. Figure 12 displays the decorated roof 



R=(P ,P 1 ,P 3 ,P 6 ,P 7 ,P 9 ,P W ,Pl2) 

of the 12— near-edge 

'0123456789 10 11 12 



E 



1010101010 1 1 
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Figure 12: A decorated roof and all its A— and V— successors 
and all its A— and V— successors 

(Po, Pl,Ps, PiiPyjPjj P9,PlO, Pn) 
(Po, Pl,P^i P&, P7, P8, -P9, -PlO) Pn) 
(Po,PuPs, Pq, P7, P9, PlO, Pll,Pl2) 

(Po, Pi,P&, Pi, P9 , Pio,Pn) 
(Po, Pi,P$, P&, P9, P10, P 12 ) 
(-Po, Pi,P^i P&, P7, P9, P12) 

A decorated roof R is a successor of a decorated roof R if R is either a 
A— or a V— successor of R. We call R an immediate successor of R if the 
difference between R and R is a minimal triangle. 

Given a finite configuration C C R 2 , a (maximal) roof-sequence is a 
sequence 

, . . . , 

of decorated roofs where Ri l is initially decorated and where Ri j+1 is a (im- 
mediate) successor of R j for j = 1, . . . ,1 — 1. The skyline of i?^ (respectively 
of is the initial (respectively final) skyline of the roof-sequence. 
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Proposition 7.1 The set of all (maximal) roof- sequences R^, . . . , with 
initial skyline S a and final skyline is in bijection with the set of all (max- 
imal) triangulations with vertices in C of the polygonal region enclosed by 
S a and Sw 

Proof. Given a roof-sequence (Ri 1 , . . . , R^ ) , the union of all associated 
skylines triangulates (maximally) the polygonal region enclosed by and 

Given a (maximal) triangulation of a polygonal region enclosed by two 
skylines, remove vertically the rightmost triangle which has only "sky" above 
it (and which is a minimal triangle for a maximal triangulation). Iterating 
this algorithm one gets an undecorated roof-sequence in reversed order. The 
decorations can now be added only in a unique way, except for the lowest 
skyline which we decorate initially. □ 
Example. The triangulation of the configuration shown in Figure 11 cor- 
responds to the roof-sequence 

(Pq, P2, P3 , -Plo); (Pqi Pit P?> , Pti P\q)i (Pq ,P2,P 7, -Pro), 
(-Po) P2, Pa, P7, Pip ), (Po, P4, P7, P10 ), (P0 P4, Pq, P7, P10), 

(Po, P4, P&, P7, P9, Pio), (Po, Pa, Pq, P7, P8,P9, Pio), 

(Po, P4, Pq, P%, P9, P10), (Po, Pa, Ps, P9, P10), (Po, Pa, P5, Pg, P9, P10) • 

Consider the rational vector space 1Z = 7Z(C) (only natural coordinates 
will in fact occur) spanned by the set of all decorated roofs of a finite con- 
figuration C C R 2 . Identifying a decorated roof with the corresponding 
basis element, we define the transfer-matrix T as the linear application 
T : 1Z — ► 1Z which sends a decorated roof to the sum of its successors. 
We introduce moreover the element G Tt(E) obtained by summing up all 
initially decorated roofs of length k with skyline d-C. Set Vo = and 

V l = C i + T(V l . 1 ) . 

We denote by W : K — > Q[s 1/2 ] the linear map which sends a decorated 
roof to if its skyline is different from d + C and wich sends a decorated roof 
R with skyline d + C to s 1 ength(i?)/2_ 

Proposition 7.2 The complete triangulation polynomial of a finite config- 
uration C = (Pq, . . . , P n ) is given by 

s f> fc / 2 W(V k ) . 

k=l 

The proof is a straightforward consequence of Proposition 17. II and of the 
behaviour of lenghts of successive roofs. 

An obvious modification of the above proposition counts triangulations 
of a polygonal (not necessarily convex) region squeezed between two skylines. 
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If one is only interested in the number of maximal triangulations of 
Conv(C), one considers the transfer matrix T which sends a decorated roof 
R to the sum of its immediate successors. The vector Cj is zero except 
for i + 1 = (t(5_C fl C) and G\ consists then of the initially decorated roof 
defined by the set d-C n C. We define Vb = and V- = C { + f (Vi-i). Set 
now W(R) = except if R is a decorated roof with underlying set <9+C n C 
where we set W(R) = 1. The number of maximal triangulations of Conv(C) 
is now given by 

E^) 

i 

(in fact, only one summand is non-zero). 

Given a near-edge E = (Po, . . . , P„), we introduce the vectors Cj € ^(E 1 ) 
and the transfer matrix T G End(7£(i<7)) as above (where we consider E as 
a configuration) and compute the vectors Vq = 0, 

Vi = Ci+T(yi-i), 1 <i<n-l . 

We introduce moreover the linear map W : K(E) — ► Q[s 1/2 ,t] defined 

by 

W(R) = s length(^)/ 2 piength(R) 

_ s length( J R)/2 ^^^-^^dength^-fe^ i length( J R)-fc 

for any decorated roof R of C. One shows now easily the following result. 

Proposition 7.3 The complete edge-polynomial of a finite near-edge E = 
(Po,.. .,P n ) is given by 

n 

E s k/2 W(V k ) . 

k=l 

7.1 A few algorithmic remarks 

Propositions 17.21 and 17.31 yield easily algorithms for computing complete 
triangulation and edge polynomials. 

These algorithms are of exponential complexity: properly organized their 
memory and time requirements are of order 2 n f(n) where f(n) is polynomial 
(of degree 1 for memory requirements, neglecting multiprecision problems, 
and of fairly small degree for time requirements) for a configuration of n + 1 
points. A nice feature of these algorithms is their simplicity: As an example, 
we run them by hand for the number of triangulations squeezed between two 
skylines of a configuration having 8 points and for the near-edge E c of the 
near-polygon depicted in Figure 2. 

For computational purposes, it is useful to have a simple representation 
for a decorated roof 

R = (Pj = Pqi Pii ) • • • ) P%d i Pid+i > ■ • • j Pim-i 5 Piw = Pn) 
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of a configuration (or near-edge) C = (Po, . . . , P n ). Such a roof R is for 
instance encoded by the integer l(R) € {0, . . . , n ■ 2 n ~ l — 1} defined by 

w— 1 

l(R) = d-2 n - 1 + ]T 2ifc_1 • 
k=l 

This yields an injection from the set of decorated roofs (of a given finite con- 
figuration (P , P n )) into a subset of Y2=\ k (k-i) = ( n + l)2 n ~ 2 natural 
integers < n ■ 2 ra_1 for n > 1. A natural integer 

n-2 

m = a ■ 2™" 1 + ^ a 2\ < q < n - 1, e» G {0, 1} 

i=0 

is of the form l(R) if and only if either a = or e Q _i = 1. The corresponding 
roof R is then given by 

R = {P , {P i+ i | < i < n - 2 and a = 1}, P n ) 

with decoration beg inning at P a . 

Let us also remark that many decorated roofs lead to dead ends when 
computing triangulation polynomials. Indeed, consider a decorated roof of 
the form 

R = {Pi , ■ ■ ■ , Pi k , ■ ■ ■ j Pi d , Pi d+1 > • • •) 

where Pi k € for k < d. If the partial skyline defined by the piece- 
wise linear path with segments consecutive points of (Pj , . . . , P{ k ) is not 
contained in d+C then there exists no roof-sequence containing R which 
stops with a roof having skyline d+(C). Indeed, R hits the "ceiling" d+(C) 
at Pi k before reaching the decorated edge. This freezes the partial skyline 
(Pj , . . . ,Pj fc ) in all successors (and iterated successors) of P. 

If one is only interested in computing the number of maximal triangula- 
tions or the maximal edge-polynomial of a near-edge, a few simplifications 
take place: 

For non-generic configurations, only decorated roofs which are full have 
to be considered where a roof is full if its sequence contains all points of C 
lying on its skyline. 

As already mentionned, the suitable transfer matrix in this case is defined 
by considering only immediate successors of a decorated roof (recall that a 
successor is immediate if it is obtained by adjoining a suitable triangle which 
is minimal, i.e. which intersects C only in its three vertices). 

Every decorated roof appears at most in one vector of Vi, . . . , V n . This 
makes memory requirements somewhat smaller. 

Let us also add that the algorithm which computes the complete trian- 
gulation polynomial of Conv(C) can be modified for computing the number 
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of triangulations (with vertices in C) for any polygonal region enclosed by 
two suitable non-crossing skylines S a , S^- Replace d±(C) by S a and S^. 

The resulting algorithm for computing triangulation polynomials or edge- 
polynomials is straightforward and best illustrated by examples. This will 
be done in the 2 following subsections. 

A few further tricks may be used: 

Computing Vi = Ci + T(Vi-±) needs only Cj and V{-\. This can be 
used to save memory by computing the contribution s^ +1 ^ 2 W(Vi-i) or 
s {i-i)/2W(Vi-i) (in the case of an edge polynomial) of V*_i to the final 
result and by erasing V^_i after determination of Vi. 

Multiprecision problems can be avoided using the Chinese Remainder 
Theorem. Compute all cofficients of decorated roofs modulo several primes 
and use them to reconstruct the result (this needs however the computation 
of an approximation of the final result which can for instance be done by a 
computation using the type of floating reals). 

Decorated roofs should always be totally ordered (using for instance the 
complete order coming from the labelling l(R) for a decorated roof R). Since 
there are at most ~ n2 n ~ 2 such roofs, the coefficient of a decorated roof 
can then be accessed using roughly 0(n log(n)) operations using a classical 
divide and conquer stategy. 

7.2 A triangulation polynomial 



We want to compute the complete triangulation polynomial of the (non- 
convex) polygonal region enclosed by the two skylines Sa^S^ defined by 
S a = (-Fb)-Pi,-FV) and Su> = (Pq, P2, P3, P5, P6, P7) of the configuration 



depicted in Figure 13. 

The following lists contain all relevant data. First we establish the list of 
all needed decorated roofs, enumerated and denoted using conventions and 
notations introduced in the previous subsection. For concision, we indicate a 
decorated roof R = (Pq, . . . , Pi d ,Pi d+1 , . . . , P7) by the (decorated) sequence 
(0 i\ ■ ■ - id id+i • • • 7) of the corresponding indices. For every roof R, we 




Figure 13: A configuration enclosed by a polygonal line 




1112 2 
13 2 12 1 
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indicate also its image T(R) obtained by applying the transfer operator T, 
except that roofs in T(R) leading to dead-ends (like (01367) for instance) 
have been omitted. For a roof R with skyline S u , we indicate s l en gth(.R)/2 
instead of T(R) which is always zero in this case. 

Data appearing in boldface can be omitted if one is only interested in 
the number of maximal triangulations. 





(07) 




Ri 


(017) 


Rn + RfiQ + i?7^+R«l + R<57 




(037) 


Rq 


Rq 


(0237) 


Rl50 + -^166 


Rs 


(047) 


Ro + R12 


Rq 


(0T47) 


Rs + -R77+Rl53 + Rl69 


R\2 


(0347) 


Ru 


i?14 


(02347) 


i?134 + R-222 + -^238 


R20 


(0357) 


R22 


R22 


(02357) 


S 2 


R36 


(0367) 


-^38 


R38 


(02367) 


-Rl82 


R52 


(03567) 


R54 


R54 


(023567) 


s 5/2 


-K69 


(013 ( ) 


/14 + lt-ji 


R71 


(0T237) 


Rq 


R73 


(0147) 


Rg + -R77 + Rl53 + Rl69 


R 77 


(0T347) 


R12 + R79 


R79 


(012347) 


Ru 


Rsi 


(0T57) 


R85 


R85 


(0T357) 


R20 + R87 


Rs7 


(0T2357 


) R22 


R97 


(0T67) 


R101 + R113 


-R101 


(01367) 


R36 + R103 


-^103 


(012367) 


R38 


R113 


(01567) 


Rin 


Rin 


(013567) 


R52 + -R119 



Rug (0123567) #54 

-R134 (0237) R150 + -R166 

R150 (02357) s 2 

R153 (01457) R 81 

R im (02367) R 182 

R im (01467) R 185 

R 182 (023567) s 5 / 2 

R 185 (014567) R 113 
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R 222 (02345 7) R 150 

#238 (023467) i?i 66 + i?254 

i?254 (0234567) i?i 82 



The vectors Cj are all except C 2 = R\. We have then V± = 0, V 2 = 
C 2 = Ri and Vi = T(Vi-i), i = 3, . . . , 9 where T is the transfer operator. 

V 2 = C 2 = i?l 

V3 = Ro + R69+-R73+R-81 + 

V4 = R4 + (R4 + R7i) + (i?8 + -R77+Rl53 + -Rl69)+R85 + (RlOl + Rll3) 

= 2R4+i?8+R7i+i?77+Rs5 + RlOl + Rll3 + R-153 + #169 

V 5 = 2R 6 +(Ro + R 1 2)+B. 6 + (R 12 + R 79 )+(Ii 20 + -Rs7) 

+ (R36 + Rl03) + Rll7 + R81 + #185 

= i?o+3R6+2-Rl2+R20 + R36 + -R79+R-81 
+R87 + Rl03 + Rll7 + #185 
Vq = i?4+3(Ri5o + Rl66) + 2i?i4+R 22 + R38 + #14+R85 
+R22 + R38 + (R52 + Rll9) + Rll3 

= R4 + 3i?i4+2R 22 + 2R38 

+R52 + R85 + -Rll3+R-119 + 3Rl50 + 3Ri66 
V7 = #6 + 3(i?i34+R 222 + i?238) + 2 ■ + 2Rig 2 

+R54 + (R20 + R87) + #117+R-54 + 3 • + 3Ri82 
= i?6+R-20 + 2R54 + Rg7 + #117 

+3i?i34 + 5Ri82 + 3R 222 + 3i?238 
V% = (Rl50+-Rl66)+R-22 + 2 • + R 22 + (i?52 + #119) 

+3(Ri50+#i66)+5 • + 3R150 + 3(i?i66 + #254) 

= 2R 22 +i?52 + i?iig + 7Ri5o + 7R\qq + 3#254 

Vg = 2 • 0+# 54 + # 54 + 7 • + 7#i 8 2 + 3#i 8 2 

= 2^54 + 10^182 

V10 = 

We have W{V 2 ) = ... = W{V 5 ) = and 



and the complete triangulation polynomial of the polygonal region between 
S a and S w is given by 



s(s 3 5s 2 + s 7/2 7s 5/2 + ^4 Qg 2 + s 9/2 ^7/2) = ^8 + ^7 + 5s 6 



7.3 An edge-polynomial 

We consider the 5— near-edge E = E c represented in Figure 14 of the near- 
polygon appearing in Figure 2. This near-edge E = (Po,...,Ps) can be 
represented by 



(W(V 6 ), . . . , W(V 9 )) = (5s 2 , 7s 5 / 2 , 9s 2 , 12s 7 / 2 ) 
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Figure 14: The near-edge E c 



In the following table we list all decorated roofs arising during the com- 
putations. The first two columns list all 23 decorated roofs needed for the 
computation, the third column yields the image T(Ri) (where T denotes the 
transfer matrix) of the decorated roof Ri described by the previous columns. 
The 6— th line for instance lists the roof (Pq, Pi, P3, P5) encoded by (0135) 
and wich we call R§ using the labelling introduced previously. An easy 
computation shows that 

T(R 5 ) = (015) + (01235) + (01345) . 

Using our numerotation of decorated roofs, we have thus 

T{R 5 ) = R 17 + P23 + R45 ■ 

Data appearing in boldface is only involved in the computation of the com- 
plete edge-polynomial and can be omitted for the computation of the max- 
imal edge-polynomial. 
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Ro (05) (O15)+(025) + (045) Ri+i? 2 + # 8 

Ri (015) (0145) R 25 

R 2 (025) (0125) + (0245) R 3 + R 26 

R 3 (0T25) (015) + (01245) R 17 + R43 

R 4 (035) (05) + (UT35) + (0235) + (0345) i? + #5 + # 6 + #28 

i? 5 (0135) (015)+(OT235) + (01345) Ri7+#23 + #45 

R 6 (0235) (025) + (DT235) + (02345) R 1S + R 7 + i? 4 6 

i? 7 (0T235) (0125) + (012345) R 35 + # 63 

Rg (045) (OT45)+(0245) Rg+#io 

R 9 (0145) 

R w (0245) (UT245) R n 

R u (UT245) (0145) R 25 

Rn (015) (0145) i?25 

R 18 (025) (0245) #26 

i? 23 (01235) (0125) + (012345) #35 + #63 

R 25 (0145) 

R 26 (0245) 

R 28 (0345) (045) i? 8 

i?35 (0125) (015) + (01245) #17 + #43 

i?43 (01245) (0145) i?25 

i?45 (01345) (0145) R 25 

ii 4 6 (02345) (0245) R 2% 

R 63 (012345) (01245) R43 



The vectors are except C 2 = R 4 . 
Vx = Ci = 

^2 = T(y 1 ) + C 2 = C , 2 = J R4 

y 3 = T(y 2 ) + c 3 = # + #5 + #6 + #28 

V 4 = T(V 3 ) + C a = (K 1 +R 2 + R 8 ) + ('R 1 7+R 2 3 + Ra5) 

+ (#18 + R7 + #46) + #8 
= Ri+i? 2 + R7 + 2i? 8 +Ri7 + i?i 8 + i?23 + #45 + #46 
V 5 = T{V 4 ) + C 5 = R 2 5 + (#3 + #26) + (#35 + #63) + 2(R 9 + J R 10 ) 

+R25 + #26 + (#35 + #63)+R25 + #26 
= i?3+2Rg + 2i?io+3R25 + 3i?26 + 2#35 + 2i?g3 

Ve = T(y 5 ) + C 6 = (R 17 + R43)+2 ■ + 2R n +3 -0 + 3-0 
+2(R 17 + # 43 ) + 2i? 43 
= 2R U + 3i?i 7 + 5i? 43 
V 7 = T(V 6 ) + C 7 = 3R 25 + 5R 25 + 2# 25 

= 10#25 
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p(E) = Y, si/2 W(Vi) 

i 

= P2S 2 + (pis 2 + 3p 3 s 3 ) + (6p 2 s 3 + 4p 4 s 4 ) 

+(13p 3 s 4 + 2p 5 s 5 ) + (3p 2 s 4 + 7p 4 s 5 ) + Wp 3 s 5 

= (10p3 + 7p 4 + 2p 5 )^ 5 + (3p 2 + I3p 3 + 4p 4 )s 4 
+(6p 2 + 3p 3 )s 3 + (pi + p 2 )s 2 
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